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ABSTRACT 

In this paper, we present the results of CARMA continuum and spectral line observations of the 
dense core Per-Bolo 45. Although this core has previously been classified as starless, we find evidence 
for an outflow and conclude that Per-Bolo 45 is actually an embedded, low-luminosity protostar. We 
discuss the impact of newly discovered, low-luminosity, embedded objects in the Perseus molecular 
cloud on starless core and protostar lifetimes. We estimate that the starless core lifetime has been 
overestimated by 4-18% and the Class 0/1 protostellar lifetime has been underestimated by 5-20%. 
Given the relatively large systematic uncertainties involved in these calculations, variations on the 
order of 10% do not significantly change either core lifetimes or the expected protostellar luminosity 
function. Finally, we suggest that high resolution (sub)millimeter surveys of known cores lacking 
near-infrared and mid-infrared emission are necessary to make an accurate census of starless cores. 
Subject headings: Astrochemistry; Stars: formation; ISM: jets and outflows 



1. INTRODUCTION 

Starless cores are the eventual birthplaces of stars, in- 
dividually or possibly in small multiples. They represent 
the transition between a diffuse molecular cloud and the 
next generation of stars to form therein. Dense cores are 
often identified by their dust continuum emission, and 
their status as either protostellar or starless can be classi- 
fied by the respective presence or absence of an embedded 
infrared object or molecular outflow. For example, cores 
in nearby molecular clouds id entified by their 1. 1 mm 
emission have been classified bv lEnoch et al.l (|2008l ). and 
cores ide ntified by their 850 nm . emission have be e n clas - 
sifie d bvlHatchell et al.l (l2007al) , Uorgensen et~aH (|2007f ) , 
and lSadavoy et al.l (|2010aT l. Given high amounts of ex- 
tinction from dust, near-infrared observations can some- 
times miss deeply embedded protostars, leading to a mis- 
classification of dense cores. Reviews of the evolution o f 
low- mass cores can be found inlDi Francesco et al.l ([2007D 
and lWard-Thompson et~aH (|2007D . 

The first "star less" core ob s erved with the Spitzer 
Space Telescope (jWerner et al.l I2004D bv the Cores to 
Disks legacy program (c2d: lEvans et al.1 120031 ). L1014, 
turned out to host an embedded infrared point source 
(|Young et al.l |2004). Other supposedly starless cores, 
such as L1448 and L1521F, were also found to harbour 
Very Low Luminosity Objects (VeL LOs) that are ei- 
ther protostars or proto-br own dwarfs (jKauffmann et all 
12003 iBourke et all f2006h . VeLLOs are defined to 
be objects embedded in den se cores that have lu- 
minosities Li n t < 0- 1 L© (|Kauffrnann et al.l 120051 : 
iDi Francesco et aLI 120071) . A search for VeLLOs in the 
c2d survey identified 50 such objects, and 15-25% of the 
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"sta rless" cores in the sur vey turned out to be misclassi- 
fied (|Dunham et"al"1 I2008D . Furthermore , iDunham et al.l 
(2008) postulated that there is likely to be a popula- 
tion of embedded objects even fainter than the sensitiv- 
ity limit of the c2d survey, leading to the conclusion that 
there is still a significant number of hidden VeLLOs. 

Despite the progress made with Spitzer, deeply em- 
bedded protostars will still be missed by near-infrared 
large-field surveys. For example, the well-studied pro- 
tostar IRAS 16293-2422 has not been detected at wave- 
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lengths shorter than 12 /mi (jJorgensen et al.ll2Qi 
cent observations of dense cores with (sub)millimcter in- 
terferometers and deeper mid-infrared observations to- 
wards known cores have found three new VeLLOs or 
candidate first hydrostatic cores in the Perseus molec- 
ular cloud. A first hydrostatic core is a short-lived 
stage that occurs after collapse has begun, between 
the prestellar core and Class protostar phases, in 
which there is hydrostatic balance between the ther- 
mal pressure of molecular hydrogen and gravity. Since 
first hydrostatic cores qualify as VeLLOs, in this pa- 
per we will group these two classes of objects together 
unless there is observational evidence t hat a particu- 
lar Ve LLO is a first hydrostatic core. IDunham et alJ 
(|2011f ) have estimated that there should be between 1- 
6 first hydrostatic cores in the Perseus molecular cloud, 
based on the expecte d lifetime of first hydrostatic cores 
(5 x 10 2 - 5 x 10 4 years: iBoss fc Yorkejl99a lOmukaill2007l 
iSaigo et all 120081 iTomida et alJl2010f >. an assumed life- 
time of the embed ded protostellar phase (5.4 x 10 5 years; 
lEvans et al.ll2009l ). and the number of embedded protot- 
stars in Perseus (66: lEnoch et al. 2009). Three candidate 
first hydrostatic cores have been identified so far in the 
Perseus molecular cloud. The first, L1448 IRS2E, previ- 
ously believed to be starless based on non-detections with 
Spitzer from 3.6 //m to 70 /im, was rec ently found to hos t 
a protostar or first hydrostatic core (|Chen et alJ 12010). 
based on the detection of a collimated CO (2 -1) outflow. 
The second, Per-Bolo 58 (jEnoch et al.l 120061 ). previously 
believed to be s tarless based on Sp itzer observations from 
the c2d survey (E vans eTalT 2003L was found to host an 
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embedded protostar or first hydrostatic core based on a 
deep 70 iim map and 3 mm continuum CARM A observa- 
tions ([Enoch, et alj|2010t iSchnee et al.ll2010D . The pres- 
ence of an embedded source in Per-Bolo 58 was then con- 
firmed with SMA observations sho wing a 1.3 mm contin - 
uum source and CO (2-1) outflow ([Dunham et al.|[20TTI ). 
Third, L1451-mm was also considered starless based 
on non-detections in the Spitzer near-infrared and mid- 
infrared images taken as part of the c2d survey. Follow- 
up observations of L1451-mm with the SMA, CARMA, 
and VLA found compact 3 mm and 1 mm continuum 
emission coinciding with broadened spectral lines (N2H + , 
NH 2 D, a nd NH 3 ) and a slow , poorly collimated CO (2-1) 
outflow (jPineda et al.ll201lD . The discovery that L1448 
IRS2E, Per-Bolo 58, and L1451-mm are not starless sug- 
gests that other "starless" cores have been misidentificd 
and instead harbour deeply embedded VeLLOs that were 
not detected in wide-area near-infrared and mid-infrared 
surveys. 

In this pap er, we present evid ence that the dense core 
Per-Bolo 45 (jEnoch et al.l 120061 despite non-detections 
in the 3.6-70 /jm range, is also protostellar. Per-Bolo 
45 is located in the Perseus molecular cloud, near the 
ce nter of the cluster N GC 1333 and between IRAS 
1 (I Jennings et all ll~987l ) to the northeast and SVS 13 
(|Strom et alJll976T ) to the southwest. In ^we describe 
new CARMA spectral line observations of Per-Bolo 45. 
The evidence that Per-Bolo 45 is protostellar is presented 
in SJ3] Some implications of the discoveries of embedded 
sources within Per-Bolo 45, Per-Bolo 58, L1448 IRS2E, 
and L1451-mm are discussed in 21 and our results are 
summarized in [JSJ 

2. OBSERVATIONS 

Here we present spectral line maps of the dense core 
Per-Bolo 45 made with the Combined Array for Research 
in Millimeter- wave Astronomy (CARMA). We made new 
observations of the 3 mm transitions of HCN, NH2D, SiO, 
HCO+, HNC, N 2 H+, and C 34 S at -5" resolution, as 
described in ij2.ll We discuss these observations in the 
co ntext of the 3 mm con tinuum data previously published 
bv ISchnee etaLl (|2010l) and described briefly in fgH 

2.1. New CARMA Observations 

Spectral line observations in the 3 mm window were 
obtained in July 2010 with CARMA, a 15 element inter- 
ferometer consisting of nine 6.1 meter antennas and six 
10.4 meter antennas. The CARMA correlator records 
signals in eight separate bands, each with an upper and 
lower sideband. We configured one band for maximum 
bandwidth (495 MHz with 95 channels per band) to ob- 
serve continuum emission, providing a total continuum 
bandwidth of approximately 1 GHz. The other seven 
bands were configured to observe spectral lines and have 
383 channels covering 7.8 MHz per band, providing na- 
tive spectral resolution of 20.4kHz (0.067 km s -1 ), which 
we smoothed to O.lkms -1 before analysis. The obser- 
vations were centered around 90.7 GHz, and range from 
85.9 GHz to 97.9 GHz. The half-power beam width of 
the 10.4m antennas is 74" at the observed frequencies. 
Seven-point mosaics were made toward the dense core 
Per-Bolo 45 (with the phase center located at J2000 
3:29:07.7 +31:17:16.8, and rest velocity 8.483 kms -1 ) in 



the D and E-array configurations, with baselines that 
range from 11m to 150 m. The largest angular scale 
to which the observations are sensitive is approximately 
30". A summary of the CARMA observations is pre- 
sented in Table [II 

The observing sequence for the CARMA observations 
was to integrate on a phase calibrator (0336+323) for 
3 minutes and Per-Bolo 45 for 21 minutes. The total 
time spent on this project was about 8 hours in each 
of the array configurations. In each set of observations, 
3C84 was observed for passband calibration, and obser- 
vations of Uranus were used for absolute flux calibration. 
Based on the repeatability of the quasar fluxes, the esti- 
mated random uncertainty in the measured source fluxes 
is a ~ 5%, and the systematic uncertainty is approxi- 
mately 20%. Radio pointing was done at the beginning of 
each track and pointing constants were updated at least 
every two hours t hereafter, using eith er radio or optical 
pointing routines ([Corder et al.l l2010h . Calibration and 
imaging were done using the MIRlAD data reduction 
package (|Sault et al.lll995fl . 

Due to instrumental problems, the first correlator band 
was lost, so we have no data on the HCN (1-0) line. We 
detected no significant emission of the C 34 S (2-1) line, 
and it will not be discussed further in this paper. 

2.2. Previous Observations 

Per-Bolo 45 was included in a CARMA 3 mm contin- 
uum surve y of 11 starless cor es in the Perseus molec- 
ular cloud (|Schnee et al.l l2010f ). It is one of two cores 
that were detected, and has an integrated 3 mm flux den- 
sity of 11+0.5 mJy, corresponding to a mass of 0.8 M© 
after making reaso nable assumptions f or the core and 
dust properties. See ISchnee et al.l (j2010l ) for more details 
about the 3 mm continuum observations and the derived 
core properties. 

Per-Bolo 4 5 wa s cla ssified as a starless core by 
lEnoch et all ([20061 120081 ). If a protostar was found 
within one full-width half-maximum (FWHM) of the cen- 
ter of the 1.1mm Bolocam core (i.e., 31", or 0.04 pc), 
then it was classified as protostellar, otherwise the 
Bolocam core was classified as starless. Protostars in 
lEnoch et al.l ([2008) were identified by the shape of their 
near-infrared to far-infrared spectral energy distributions 
(SEDs), a minimum flux at 24 /im and the presence 
of 70 /im point sources not classified as galaxy candi- 
dates. Per-Bolo 45 was iden tified as a Class proto- 
star by lHatchell et al.1 (|2007al ) based on the colors of a 
nearby, but (we believe) unrelated, near-infrared source. 
Th ere is a possib l e outflo w from Per-Bolo 45 reported 
by lHatchell et al.l (|2007bD . but due to confusion from 
nearby p rotostars this id e ntifica tion is given as tenta- 
tive. In iJorgensen et al.1 (|2007h . cores were identified 
as protostellar based on either the presence of Spitzer 
detections and colors within a SCUBA core, or by the 
concentratio n of the light pro f ile of t he 850 /im SCUBA 
flux Neither |j0rgensen et all (|2007D . Kirk et al.l (|2007t) . 
nor iDunham et al.l (|2008f ) found evidence for an embed- 
ded source i n Per-Bolo 45, and it is classified as a star- 
less core by iSadavov et ail (|2010af ). In Fig. [U we show 
Spitzer near-infrared and mid-infrared images of Per- 
Bolo 45 with CARMA 3 mm continuum overlays, and 
it is clear that there is no embedded point source de- 
tected by Spitzer associated within Per-Bolo 45. The 
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area shown in Figures [T] - [5] is approximately the re- 
gion in which the gain of the mosaic is greater than 0.5. 
Therefore, the true classification of Per-Bolo 45 is not 
clear from the literature, but most studies have labeled 
it "starless." 

3. ANALYSIS 

We study the kinematics and morphology of the dense 
core Per-Bolo 45 by fitting the line profiles of NH2D (lu 
- 1 1), SiO (2-1), HCO+ (1-0), HNC ( 1-0), andN?H+ (1 - 
0) using the MPFIT package in IDL (|Markwardt)l2009[ ). 
Single Gaussians were fit to the SiO, HCO + , and HNC 
spectra, and multiple Gaussians were fit to the 7 com- 
ponents of the N 2 H + spectra and 6 components of the 
NH 2 D spectra. Fits were only made to those profiles that 
had three independent velocity channels with a signal-to- 
noise ratio greater than 3. 

3.1. Molecular Line Emission 

Starless cores are kinematically quiescent objects, with 
line widths broadened approximately equally by ther- 
mal and non-thermal motions. The thermal line widths 
(FWHM) for the molecular lines presented in this pa- 
per a re ~0.13 km s -1 , as suming a typical temperature of 
10 K (|Schnee et al.ll2009h . The average non-thermal line 
width for NH3 in starless cores in the Perseus molecular 
cloud is greater than the thermal line width of NH3 by 
a factor of 1.5 (|Foster et alj l2009h . Broad line widths 
(FWHM - 0.6 km s" 1 ) in N 2 H+ and NH 2 D towards the 
3 mm continuum peak of L1451-mm led I Pineda et al.1 
(|201 If ) to hypothesize the presence of an embedded pro- 
tostar in a core that had previously been considered 
starless, which was confirmed by detection of a CO 
(2-1) outflow. Similarly, velocity gradients in starless 
cores are small, with a typical value of 1 kms _1 pc _1 
(|Goodman et al.l [1993). We would expect a total ve- 
locity gradient across Per-Bolo 45, which has a diame- 
ter no greater than 0.2 pc, to be ~0.2kms _1 . Further- 
more, we would expect the velocity field to be continu- 
ous, a s seen in interfero metric ma ps of the starless core s 
L183 (TKirket al.l [2009h. L694-2 (jwilliams et "all [20061) . 
and L1544 ([Williams et al.ll2006D . The kinematics of Per- 
Bolo 45, as traced by N 2 H+ N H 2 D, H CO+, and HNC, 
are described in Section l3.1.1l and l3.1.2l Maps of the full- 
width half-maximum (FWHM) of N 2 H+, NH 2 D, HCO+, 
and HNC are shown in Figure [U and the LSR velocity 
maps for N 2 H+, NH 2 D, HCO+, and HNC are shown in 
Figure [3] Example spectra demonstrating the variety of 
line profiles in Per-Bolo 45 are shown in Figure [4] 

In the dense, cold, and quiescent interstellar medium, 
SiO is frozen onto the surfaces of dust grains except 
in regions shocked by outflows t hat return SiO to the 
gas p hase and excite its emission ([Martin-Pintado et al.l 
1992). The presence of SiO emission is therefore a good 
indicator of one or more nearby protostars. SiO emis- 
sion can be grouped into two classes, one that has broad 
lines (>10kms _1 ) and is Doppler-shifted by velocities 
^lOkms^ 1 from the ambient cloud, and a second com- 
ponent with narrow lines (~1 kms -1 ) found at velocities 
very close to the cloud ma terial (e.g.. lLefioch et al1l!998l : 
Uimenez-Serraet al.|[200l ). The SiO maps of Per-Bolo 
45 are discussed in Section [3.1.31 and presented in Figure 
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An 850 fxm SCUBA map of the dense core Per-Bolo 
45 was produced as part of SCUBA Legacy Catalogues 
(|Di Francesco et aTll2008[ ). and is shown in Figure[6]along 
with the 3 mm continuum emission and integrated inten- 
sity maps of N 2 H+, NH 2 D, HCO+, and HNC. Whereas 
the 3 mm continuum emission detected with CARMA has 
only one peak, the integrated intensity maps of N 2 H + , 
NH 2 D, HCO+, and HNC all have multiple peaks. The 
morphology of emission towards Per-Bolo 45 varies sig- 
nificantly between molecular lines, and this will be dis- 
cussed below and presented in Figure [6] 

3.1.1. N 2 H+ and NH 2 D 

The line width of N 2 H + in Per-Bolo 45 has a bimodal 
distribution, with narrow lines (~0.3kms _1 ) at the posi- 
tion of the dust peak and a small region of much broader 
(0.5-1.2 kms -1 ) line width a few arcseconds to the south- 
east of the dust peak. The line width of NH 2 D is qual- 
itatively similar to that of N 2 H + , with a narrow line 
width at the dust peak and much broader width a few 
arcseconds to the southeast. We argue that the broad 
line widths (up to 1.2 kms -1 ) seen in N 2 H + and NH 2 D 
towards Per-Bolo 45 point towards the presence of a pro- 
tostar, which must be low luminosity to have not been 
detected in mid- infrared Spitzer maps (see § 13. 2p . The 
velocity fields of N 2 H+ and NH 2 D (as well as HCO+ and 
HNC) are shown in Fig.[3j The range of velocities seen 
in N 2 H + and NH 2 D is fairly narrow and similar to the 
systemic velocity of the core as measured in the single 
dish maps. 

Both N 2 H + and NH 2 D are coincident with the 850 nm 
dust emission, but extend about 30" to the southeast be- 
yond the extent of the 3 mm continuum. In both lines the 
secondary peaks have narrow line widths typical of quies- 
cent cores, so we suggest that this material is part of the 
envelope of Per-Bolo 45 that is not strongly affected by 
any outflow. The secondary peaks of N 2 H + and NH 2 D 
are coincident with the elongation of the 850 /im dust 
emission seen in the SCUBA map, although they are not 
accompanied by a similar feature in the 3 mm continuum 
map. Discrepancies on small spatial scales between the 
dust emission and dense gas tracers have been seen be- 
fore, for in stance in the dense core Oph B dFriesen et all 
120091 120101) and in Oph A-N6 (|Pon et al.ll2009D . Further^ 
more, we show in this paper that the CARMA continuum 
maps are not sensitive to the emission from starless ma- 
terial. 

3.1.2. HCO+ and HNC 

The HCO + integrated intensity map (see Fig. [6]) shows 
two peaks, one immediately to the southeast of the dust 
peak and the other 40" further to the southeast. The in- 
tegrated intensity map of HNC shows a peak immediately 
to the southeast of the dust peak and also has three more 
peaks to the east and southeast of the 3 mm continuum 
emission. Two of the HNC emission peaks to the east 
and southeast have no counterparts in the dust emission 
or with the other features in molecular line maps. Both 
HCO + and HNC show a broad range of velocities and 
both have line widths many times greater than the ther- 
mal line width, as shown in Figs. [5] and G3 In the case of 
HCO + , both redshifted and blueshifted emission are seen 
with a spread of ±1 kms -1 from the systemic velocity of 
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the core, and no emission is seen at the systemic veloc- 
ity of Per-Bolo 45. The HNC emission adjacent to the 
3 mm continuum peak is at the systemic velocity of the 
core and has broad line widths (up to 2 kms -1 ), while 
the other HNC peaks are blue-shifted and have narrow 
line widths. 

Although the interpretation of the emission is com- 
plicated by dissimilarities between the morphologies of 
the emission, we argue that the broad line widths (up 
to 2 kms -1 ) and velocities Doppler-shifted from the sys- 
temic velocity of the core (by ±lkms _1 in the case of 
HCO + ) are likely caused by an outflow driven by Per- 
Bolo 45. Neither kinematic feature would be expected 
in a starless core, as explained in Section 13.11 The re- 
gions of narrow line width seen in HNC and HCO + (as 
well as N 2 H + and NH2D) likely trace the more quies- 
cent gas in a dense envelope. HCO + has been seen to 
trace outflows in the cluster NGC 1333 fe.g. JWalsh et al.l 
2007), of which Per-Bolo 45 is a member. Given that 
HCN has been also see n to trace outflows in NGC 1333 
(|j0rgensen et al.ll2004f) . it is reasonable to suggest that 
the broad line width HNC in Per-Bolo 45 traces an out- 
flow as well. HNC has been seen to trace dense ma- 
terial around a protostar as well as material associ ated 
with the protostellar outflow (|Arce fc Sargentl I2004D . so 
the appearance of both quiescent and turbulent HNC in 
Per-Bolo 45 is not unexpected if it harbors a protostar. 

3.1.3. SiO 

Per-Bolo 45 exhibits SiO emission, with ~lkms _1 
line widths (see Fig. [4]) and velocities within ~lkms _1 
of the ambient medium (see Fig. [5]). The systemic ve- 
locity of Per-Bolo 45 i s -8.5 km s" 1 (jKirk et al.l 120071: 
Rosolows kv et al J 12008). The systemic velocity of the 
ISM around Per-Bolo 45, measured in an FCRAO map o f 
CS (2-1) by the COMPLETE survey (jRidge et al]l200fih . 
is 8.6km s -1 , with a 1.7km s line width. The veloc- 
ity range of the SiO (2-1) emission towards Per-Bolo 45 
is therefore entirely consistent with the range of veloc- 
ities of the ambient material. The morphology of the 
SiO emission is perhaps suggestive, being extended with 
the long axis of the emission pointing towards the 3 mm 
continuum emission. 

We interpret the detection of SiO (2-1) emission im- 
mediately adjacent to Per-Bolo 45 as further evidence 
for the presence of an embedded protostar. An outflow 
launched by a protostar in Per-Bolo 45 would liberate 
SiO from the dust in the ISM, creating the conditions 
required for the observed emission. Since we only detect 
the narrow component of SiO emission that traces the 
ambient cloud material, the velocity map shown in Fig. [5] 
is not sufficient to say for certain what portion of the 
outflow is exciting the emission. The velocity structure 
seen in Fig.[S] is likely that of the ISM around Per-Bolo 
45, and not that of the outflow. The velocity range cov- 
ered by the CARMA observations (-5 kms -1 < VLSR < 
21 km s -1 ) was not wide enough to find the high- velocity, 
broad component of the SiO emission more closely asso- 
ciated with the outflow, if present. 

3.1.4. The case for the protostellar nature of Per-Bolo 45 

We have shown that Per-Bolo 45 exhibits several traits 
of protostellar cores, despite the non-detection of mid- 



infrared emission in IRAC and MIPS Spitzer maps (see 
Fig.HJ). The line widths seen in N 2 H+, NH 2 D, HCO+, 
and HNC are several times larger than the thermal line 
width, behaviour not found in starless cores but which 
could be explained by the presence of an embedded pro- 
tostar. Broad line widths of N 2 H+ and NH 2 D in L1451- 
mm led to the discovery of a protostar in that core 
even though it had pr eviously been identified as star- 
less (jPineda et al.ll20ll . The presence of redshifted and 
blueshiftcd HCO + emission seen about 1 kms -1 from the 
systemic velocity of the core would also not be expected 
in a starless core, but could be explained by an outflow 
from an embedded protostar. SiO emission located adja- 
cent to Per-Bolo 45 can be explained by the interaction 
between the ambient material in the molecular cloud and 
an outflow launched by the core. Finally, the 3 mm con- 
tinuum emission detected towards Per-Bolo 45 is also 
suggestive of a possible embedded source, given that 
3 mm continuum emission was detected with CARMA to- 
wards Per-Bolo 58 and L1541-mm, two supposedly "star- 
less " cores that were subsequently found to be protostel- 
lar (lEnoch et alj|201tt [Dunham et al.ll2011t|Pmeda et al.l 
120111) . Non e of the othe r starl ess cores surveyed with 
CARMA bv lSchnee etaLl pOlOh exhibited 3 mm contin- 
uum emission. 

Given that Per-Bolo 45 is in the NGC 1333 cluster, it 
is important to check that the outflow does not originate 
from another nearby protostar. The elongation of the 
SiO emission points back towards the HH 7-11 group as- 
sociated with SVS 13, and even further perhaps towards 
IRAS 4, but each of these cores have well-known outflows 
going in directions other than towards Per-Bolo 45. For 
example, a CO (3-2) map of NGC 1333 shows no out- 
flow directed towar ds Per-Bolo 45 from any other object 
(jCurtis et al.l 12010V Instead, Per-Bolo 45 seems to be 
adjacent to a compact red lobe emanating from IRAS 7 
to the northeast and a very extended red lobe emanat- 
ing from IRAS 2 in the southwest. These redshifted lobes 
have velocity ranges of 12 kms -1 to 18 kms -1 and their 
blueshifted lobes range from -5 kms to 3 kms -1 , both 
outside the range of velocities associated with the red- 
shifted and blueshifted HCO + emission from Per-Bolo 
45. It would also be difficult to explain how an outflow 
from elsewhere in the NGC 1333 cluster could interact 
with the dense material in Per-Bolo 45 traced by N 2 H + 
and NH 2 D in such a way as to increase turbulence im- 
mediately to the south-east of the 3 mm continuum peak 
without also affecting the area with nearly thermal line 
widths around it. We conclude that the molecular line 
observations presented in this paper are best explained 
by the presence of a protostar embedded within Per-Bolo 
45. 

3.2. Embedded Source 

A relationship between the internal luminosity of a 
young stellar obj ect (YSO) and i ts 70 /mi flux is given 
by Equation 2 in iDunham et al.l (|2008T > . Per-Bolo 45 is 
not detected at 70 /xm, but we can use a 3 a upper limit 
to its 70 /im flux to estimate an upper limit to the inter- 
nal luminosity of the embedded source. We find that the 
internal luminosity of Per-Bolo 45 is less t han 10~ 2 Lp), 
simila r to the upper limit for L1451-mm bv lPineda et al.l 
(|2011[) and lower by a factor of about 10 than the lumi- 
nosities of the embedded protostars observed by Spitzer 
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with p ublished models shown in Table 1 of lDunham et all 
(2008). The upper limit to the luminosity of the embed- 
ded source in Per-Bolo 45 is consistent with the sensi- 
tivity limit of the c2d survey. Assuming that Per-Bolo 
45 is protostellar, the lack of observed 70 /jm flux would 
make this source a VeL LO by the definition given in 
iDi Francesco et ail ([2007D . A protostellar source with 
a disk viewed edge-on might remain invisible at the ob- 
served sensitivity limits, so it is possible that the inter- 
nal luminosity of Per-Bolo 45 is higher than that of a 
VeLLO and it appears dim because of the viewing angle. 
Given its low 70 /jm luminosity, the source embedded in 
Per-Bolo 45 is also a plausible first hydrostatic core can- 
didate. 

4. DISCUSSION 

In this paper we identify Per-Bolo 45 as a core with an 
embedde d source, changing its classification from star- 
less core (jEnoch et alj 12008) to VeLLO. Other cores m 
the Perseus molecular cloud, previously identified as star- 
less due to their lack of near- and mid-infrare d emission 
(e.g. lEnoch et al] 120081 : ISadavov et all l2010aft but sub- 
sequently found to show evidence for protostellar activ- 
ity, such as launching molecular outflows, incl ude Per- 
Bolo 58 (lEnoch et al.ll2010l : lDunham et allboill) . L1448 
IRS2 E (|Chen et al.ll2010D . and L1451-mm ( fPineda et al.l 

EoTD). 

Out of the 11 "st arless" cores in Perseus surveyed by 
iSchnee e"t~ al. (2010) with interferometric observations of 
their 3 mm continuum, the only two cores detected (Per- 
Bolo 45 and Per-Bolo 58) have since been shown to har- 
bour embedded objects. This fraction of misidentificd 
starless cores in Perseus, 2/11 or 18% , is probably an up - 
per limit given that the cores in the ISchnee et all (|2010h 
sample were chosen to have hi gh surface brightnes ses in 
the 1.1mm continuum map of lEnoch et al. (2006), and 
the presence of a protostar would increase the peak flux 
of a dense core. 

A lower limit to the number of cores in Perseus 
misidentified as starless can be derived from the number 
of previously identified starless cores and the number of 
newly identified protostars and VeLLOs. This estimate 
provides a lower limit because not every core identified 
in wide-field dust continuum maps of Perseus has been 
followed up with (sub)millimeter interfer ometric molec- 
ular lin e observations to find outflows. Sa davoy et all 
(|2010al ) find 97 starless cores and 46 protostellar cores 
in Perseus, using SCUBA 850 /im maps to identify the 
cores and Spitzer maps between 3.6 and 70 /xm to de- 
termine the starless o r pro tostellar status of each core. 
!J0rgensen et al.l (|2007l ) and lEnoch et all (|2008l ) find sim- 
ilar numbers of starless cores and protostars using simi- 
lar techniques. Although the census of protostellar cores 
and starless cores detected by Herschel has not yet been 
released for the Perseus molecular cloud, a 2:1 ratio of 
prestellar co res to protostell a r core s in t he Aquila rift was 
report ed by IKonvves et al.l (|2010() and iBontemps et all 
(2010), in good agreement with the pre- Herschel studies 
of Perseus described above. Considering the four newly- 
identified VeLLOs in Perseus, a lower limit of 4/97 "star- 
less" cores, or 4%, have been misclassified. Our upper 
limit of 18% would imply that 17 cores have been mis- 
classified as starless. We note that it is likely that not all 
of the dense cores in Perseus have been identified, given 



that the published (sub) millimeter surveys used to iden- 
tify such cores have had limited sensitivity to faint and 
extended structures. 

The lifetime of starless cores can be estimated from the 
relative number of starless cores and protostellar cores, 
given a reasonable estimate of the protostellar lifetime 
and assuming t hat the rate of sta r formation is con- 
stant over time. lEnoch et all (|2008l) found that a mean 
lifetime for starless cores in nearby molecular clouds is 
0.5±0.3Myr, so the <20% change in the number of star- 
less cores in Perseus resulting from the discovery of em- 
bedded sources in previously identified "starless" cores 
is roughly comparable to the uncertainty coming from 
other systematics in the estimate of starless core life- 
times. 

Similarly, increasing the number of protostars in c reases 
the resultant protostellar lifetime. lEvans et all (2009) 
found 87 Class 0/1 protostars in the Perseus molecular 
cloud, so adding 4-17 more sources into the Class 0/1 cat- 
egory would increase the Class 0/1 lifetime by ^5-20%. 
This increase in the Class 0/1 lifetime would not signif- 
icantly affect modelling of the expected accretion lumi- 
nosity averaged over the lifetime of a protostellar core, 
given the significant uncertainties on both the observa- 
tional a nd theoretical sides of the problem, as summa- 
rized in lOffner fc McKeel (|20lTI ) . 

Many studies o f nearby molecular clouds (e.g. , 
i Motte et all 119981: iNutter fc Ward-Thompsonl 1 2007 : 
Ward-Thompson et all [2007E lEnoch et alf 12008 : 
IKonvves et al.l I2010D have shown that the starless 
core mass function (CMF) has a similar shape to 
the stellar initial mass function (IMF). These studies 
depend on the ability to disti nguish starless den s e core s 
from protostellar d ense cores. lHatchell fc Fullerl (I2008T) . 
lEnoch et all (|2008l ) and IDi Francesco et al.l (|2010f ) each 
found that their most massive cores tended to be 
protostellar, which suggests that either prestellar cores 
have an upper mass limit, consistent with gravitational 
stability inequalities (e.g., Jeans analysis), or that the 
more massive prestellar cores have ve ry short lifetimes. 
Furthermore, ISadavov et al.l (|2010bl ) found that the 
most massive "starless" cores in Orion and Perseus had 
ambiguous infrared emission towards them and thus, 
changed their classification to "undetermined." If the 
most massive "starless" cores tend to be misclassified 
VeLLOs, then the CMF will lack the massive core 
tail found in the IMF. This difference will steepen the 
starless CMF slope with respect to the IMF slope, i.e., 
the two mass functions will have different shapes. 

Three of the VeLLOs in Perseus (Per-Bo lo 45, Per-Bolo 
58, and L1448 IRS2e) were included in the lSadavov et all 
(|2010bf ) investigation into the gravitational stability of 
starless cores in Perseus. None of these sources were 
found to be particularly unstable. Instead, they reside, 
along with the majority of starless cores, in a mass range 
of 1-3 Jeans masses (see ISadavov et al1l2010bl for details 
on how the stability analysis was performed). The only 
possible hint that these cores are unique is that all three 
have relatively small sizes (compared with other cores 
at similar level of Jeans instability). That these cores 
are not clearly distinguishable by a Jeans analysis fur- 
ther suggests that deep mid-infrared and interferometric 
observations will be required to uncover the majority of 
VeLLOs. Given that only a few "starless" cores have 
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been proven to be misclassified, and that these cores do 
not stand out as being especially massive, the shape of 
the CMF, as discussed above, has not been biased by the 
presence of low-luminosity protostellar objects. 

In a 3 mm contin uum survey of 11 "starless" cores, 
iSchnee et all (|2010t l found that only two (Per-Bolo 45 
and Per-Bolo 58) were detected and concluded that the 
density distribution in starless cores must be smooth 
and not strongly peaked. The su bsequent discovery 
that both cores are a ctually VeLLOs ([Enoch et al.ll2010t 
iDunham et al.ll2011l this paper) strengthens the claim 
that starless cores in Perseus have smooth and shallow 
density profiles. A shallow inner density profile in star- 
less c ores has been previously reported by s everal groups, 
(e.g iWard-Thompson et all 11994 119991: iShirlev et ail 
2000) . Had the starless cores in our Perseus sample been 
in the process of fragmentation, we would have been able 
to detect the resultant fragments in the dust emission 
maps. 

Although Spitzer surveys were more sensitive than 
those conducted with previous instruments and found 
embedded objects in cores that had been identified 
as starless (lYoung et al.l 120041 : iKauffmann et al.1 120051 : 
iBourke et al.ll2006l ) . it is also true that deeply embedded 
protostars can be miss ed at wavelengths les s than 12 iim 
(|j0rgensen et al.ll2008D by surveys li k e c2djEvans et al.l 
12003ft . In addition. IDunham et al.l (|2008| ) have found 
that there is likely to be a population of VeLLOs too 
faint to have been detected in current surveys. New Her- 
schel observations will be able to find fainter protostars 
than Spitzer was able to detect, and in the Aquila region 
seven of the ~50 Class protostars detected by H erschel 
were missed by Spitz er (|Bontemps et al.l 12010ft . Still, 
iBontemps et al.l (|201Q[ ) reported that compact sources in 
Aquila were identified down to the Herschel 70 /im de- 
tection limit, implying that a population of even fainter 
sources have yet to be discovered. Given that low- 
luminosity embedded objects can be found through high- 
resolution observations of their dust continuum emission 
and molecular outflows, we suggest that interferometric 
(sub)millimeter observations are a promising method for 
determining what fraction of "starless" cores have been 
misclassified. 

5. SUMMARY 

In this paper, we report on CARMA maps of the 3 mm 
continuum and 3 mm window spectral lines of NH2D, 
SiO, HCO+, HNC, and N 2 H+. Our main results are: 

• Despite non-detections in Spitzer maps at 3.6- 
24 /xm, Per-Bolo 45 is a protostellar core, as inferred 



from large line widths, Doppler-shifted emission, 
and the presence of SiO emission. 

There are at least four cores in Perseus, previ- 
ously identified as starless, that recent observa- 
tions have shown to contain low-lumino sity embed 



ded o bjects with m olecular outflows dChen et all 
20101: (Enoch et all 120101 IDunham et al.l 120111 : 



Pineda et alll201lF . If this result can be general- 



ized to other nearby molecular clouds, we estimate 
that the lifetime of starless cores has been overesti- 
mated by 4-18% and the lifetime of Class 0/1 pro- 
tostars has been underestimated by 5-20%. These 
changes are within the previously published uncer- 
tainties of starless and protostellar c ore lifetimes 
(|Enoch et al.ll2008HEvans et al.ll2009h . 

• Although recent infrared surveys of nearby 
molecular clouds have made great progress to- 
wards classifying star le ss and p rotostellar cores 
(e.g., lHatchell et al.l l2007al Uorgensen et al l 



2007 



2008 



20081: lEnoch et all 12001 IDunham et al.l 
Evans et al.l 120091: ISadavov et all 12010a ' 



Konvves et al.l 12010b IBontemps et all 120101 ) . we 
suggest that an improved census will require high- 
resolution (sub)millimeter observations to survey 
known cores and identify embedded low-luminosity 
objects and outflows. 
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Fig. 1. — Spitzer near-infrared and mid-infrared maps of Per-Bolo 45 (greyscalc) with CARMA 3 mm continuum contours. The Spitzer 
maps at 3.6 ^tm (t op left), 4.5 (im (top right), 8.0 fim ( bottom left) and 24 fim (bottom right) were taken as part of the c2d survey 
(Evans ct al. 2003). The 3mm continuum data come from Schncc ct al. (2010), and only the 5<r flux density contour is plotted. The (0,0) 
position is J2000 3:29:07.7 +31:17:16.8, and is the same in all figures in this paper. The units of the maps are MJy/sr. 
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Fig. 2.— The full-width half-maximum (FWHM) of N 2 H+ (top left), NH 2 D (top right), HCO+ (bottom left), and HNC (bottom right). 
The line width is given in color, as shown in the scale bar to the right of each panel, in units of km s — 1 . See JHfor an explanation of the fits 
to the spectral line profiles. Line widths are only calculated for positions with at least three channels with emission above the 3<r cutoff. 
The red contour shows the 5 a dust continuum emission, as in Fig. [T] Th e black contours show the SCUBA 850 /im map of Per-Bolo 45 at 
70% and 90% of the peak emission in the map (Di Francesco et al. 2008). The synthesized beam size of the continuum emission is shown 
in the bottom left corner of each panel, and the beam size of the spectral lines is shown in the bottom right corner of each panel. 
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Fig. 3.— The VLSR of N 2 H+ (top left), NH 2 D (top right), HCO+ (bottom left), and HNC (bottom right). The fit velocity is given in 
color, as shown in the scale bar to the right of each panel, in units of kms - 1 . See §3\ for an explanation of the fits to the spectral line 
profiles. Velocities are only calculated for positions with at least three channels with emission above the 3 <r cutoff. The red contour shows 
the 5 ex dust continuum em ission, as in Fig. [T] The black contours show the SCUBA 850 fim map of Per-Bolo 45 at 70% and 90% of the 
peak emission in the map (Di Francesco ct al. 2008). The synthesized beam size of the continuum emission is shown in the bottom left 
corner of each panel, and the beam size of the spectral lines is shown in the bottom right corner of each panel. 
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Fig. 4. — Example spectra towards Per-Bolo 45. SiO emission (top row) is presented from a relatively blueshifted position (left) and 
redshifted position (right). For HCO + , HNC, NH2D, and N2H+ (second from the top to the bottom, respectively), relatively wide (left) 
and narrow (right) line profiles are shown. The positions of the spectra are shown in the top right corner of each panel. Data are shown 
with the solid black line, and the fit profiles are shown with the dotted red line. 
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Fig. 5. — SiO (2-1) emission from Per-Bolo 45. Color indicates the VLSR of the emission and dotted black contours show the peak 
antenna temperature (1, 3, 5K). The red contour shows the 5<r dust mv cont inuum emission. The soli d black contours show the SCUBA 
850 /an map of Per-Bolo 45 at 70% and 90% of the peak emission in the map (Di Francesco et al. 2003). The synthesized beam size of the 
continuum emission is shown in re d in the bottom left corner, and the SiO beam size is shown in black in the bottom right corner. The 
VLSR of Per-Bolo 45 is 8.5 km s" 1 IIKirk et al. 2007; Rosolowskv et al. 200|). The SiO emission is consistent with gas liberated from dust 
grains and excited by an outflow launched by Per-Bolo 45. 
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Fig. 6.— The integrated intensity of N 2 H+ (top left), NH 2 D (top right), HCO+ (bottom left), and HNC (bottom right). The integrated 
intensity is given in color, as shown in the scale bar to the right of each panel, in units of Kkms . See [[3]for an explanation of the fits to 
the spectral line profiles. The red contour shows the 5 <r 3 mm dust continuum emissi on, as in Fig.JT] T he black contours show the SCUBA 
850 fim map of Per-Bolo 45 at 70% and 90% of the peak emission in the map jE)i Francesco ct al. 20Q8|). The synthesized beam size of the 
3 mm continuum emission is shown in the bottom left corner of each panel, and the beam size of the spectral lines is shown in the bottom 
right corner of each panel. 



